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Two-particle correlations of direct photon triggers with associated hadrons are obtained by isola- 
tion method in p + p collisions at y/s= 200 GeV in PHENIX at RHIC. The initial momentum of 
the away-side parton is tightly constrained, because the parton-photon pair is balanced in momen- 
tum at the leading order in perturbative quantum chromodynamics (pQCD). Therefore making 
such correlations can be used as a tool to measure the away-side parton fragmentation function. 
The direct photon associated yields in p + p collisions are compared with PYTHIA [1] and the 
effect of the kj smearing in the spectra is discussed. 
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1. Introduction 

Fragmentation functions represent the probability for a parton to fragment into a particular 
hadron carrying a certain fraction of the parton's energy Fragmentation functions incorporate the 
long distance, non-perturbative physics of the hadronization process in which the observed hadrons 
are formed from final state partons of the hard scattering process. It, like structure functions, cannot 
be calculated in perturbative QCD but can be evolved from a starting distribution at a defined energy 
scale [2]. If the fragmentation functions are combined with the cross sections for the inclusive 
production of each parton type in the given physical process, predictions can be made for the scaled 
momentum, z, spectra of final state hadrons. Small z fragmentation is significantly affected by the 
coherence (destructive interference) of soft gluons [3], while scaling violation of the fragmentation 
function at large z allows to measure a s [4]. It has been believed that the shape of the high-p^ 
trigger hadron associated xe distributions reflects the shape of the fragmentation function [5], by 
xe we mean 

PTx-PTa PTa . , ZaPTa 

xe = j = cosA0~ — — (1.1) 

Pti Pti ZtPn 

where pju PTa, Ptu PTa are the transverse momenta of the trigger and associated parton, trigger and 
associated hadrons and z a =PTa/pTa and Zt=PT\.lpTt- However, it has been demonstrated in [6] that 
the fixed momentum of the trigger particle pjt does not fix the mother parton momentum pn and 
pn varies with pj a . In this case both z a and Zt vary with pj & and this variation completely masks 
the actual shape of the fragmentation function. One of the alternative ways to explore the fragmen- 
tation function is to study the particle distributions associated to the direct photon. In the high-py 
region where the photon-production is dominated by the Compton scattering, the direct photon bal- 
ances the back-to-back quark. The associated xe distribution then approximates the fragmentation 
function of the away side jet. However, in this case there are important effects which need to be 
taken into account. One constraint comes from the fact that at low pr the direct photon production 
may be contaminated by fragmentation photons from the large number of gluonic jets produced 
[7, 8]. Fragmentation photons obviously do not balance the away-side parton. Another constraint 
comes from the fact that even when considering the leading order Compton diagram there is always 
soft QCD radiation which breaks the jet-photon momentum balance and the azimuthal collinearity. 
This non-perturbative radiation manifests itself as a non-zero value of the net parton-photon pair 
transverse momentum magnitude. Originally this soft radiation induced transverse momentum was 
attributed to the individual incoming partons and the notation of kj (parton transverse momentum) 
was introduced by Feynman, Field and Fox [9], where the colliding partons might have some initial 
transverse momentum kj with respect to the incoming hadrons. This analog of the Fermi motion 
of nucleons in a nucleus would give rise to a smearing out of the pj spectrum. There could be an 
additional smearing due to the kj associated with the fragmentation in the final state [10]. 

2. PHENIX 7T° — h and direct y—h correlation results 

The identification of direct photons is difficult due to the large number of background photons 
from hadronic decays, mostly from 71° decays. Therefore the extraction of the direct photon-hadron 
pairs per trigger yields relies on a statistical subtraction of the decay photon-hadron per trigger 
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Figure 1: PHENIX preliminary data from p + p collisions at y(s=200 GeV. Left : xe distributions of 
charged hadrons associated with trigger n° (open markers) and isolated y (closed markers) are fitted to a 
simple exponential function (Ne~ bXE ) in the available range of xe- The mean xe values extracted from the 
distributions are indicated by a solid line which goes across the xe distributions labeled as "Mean xe" ■ As 
for the dotted line, labeled as "PHENIX acceptance", the mean xe is calculated by limiting the associated 
charged hadron momenta from 1 to 4 GeV/c of PHENIX acceptance in a given trigger momentum. Right : 
The extracted exponential slope b (called "inverse xe slope" in the figure) is shown as a function of trigger 
particle momentum. 



yields from the inclusive photon-hadron pairs per trigger yields. The first result from this "sub- 
traction method" was presented in [11]. Recently PHENIX developed photon isolation cuts which 
were applied event by event in a new analysis in order to dramatically reduce the contamination of 
di-jet events with 71° decay or fragmentation photons in the y — h event sample. The measured xe 
distributions of charged hadrons associated with 71° and isolated photon triggers in p + p collisions 
at y/s = 200 GeV are presented in Fig. 1 for various isolated photon and 7T° trigger momenta. A fit 
of the xe distribution to a simple exponential (Ne~ bXB ) in the available range of xe is presented in 
left panel of the Fig. 1. The extracted exponential slope b from the fitting as a function of trigger 
pn is shown in the right side of the Fig. 1. A very interesting feature from this result is observed. 
The xe slopes are still increasing as the trigger momentum gets larger, which is very similar to the 
n° — h case. 

In order to understand this behavior of xe slopes as a function of trigger momentum, we 
compare the PHENIX preliminary data with PYTHIA, where J (kj^ = 3.0 GeV/c is used based on 
the PHENIX measurement [6]. Fig. 2 shows the xe distributions of quark and gluon fragmenting to 
71 s for various quark and gluon momenta calculated by PYTHIA (note that xe is negative because 
only near side particles associated with quarks and gluons are accounted for). The xe distributions 
of 7r°'s resulting from quark and gluon fragmentation are parameterized by means of a fit function 

^«jt i «(l-jc B )/>, (2.1) 
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(a) PYTHIA quark fragmentation function 



(b) PYTHIA gluon fragmentation function 



Figure 2: xe distributions of n s resulting from quark and gluon fragmentation in various quark/gluon 
momenta, pj t in PYTHIA (p + p at y/s=200 GeV). The results of the KKP parametrization [2] fitting to 
the distributions are shown as solid (for u-quarks) and dashed (for gluons) lines. The momenta in the figure 
label are given in GeV/c. 

i.e. we use in fitting the same functional form as in KKP parametrization [2] for fragmentation 
functions. The results of the fitting are shown as solid (for u-quarks) and dashed (for gluons) lines. 
Negative Logarithmic Derivative (NLD) of xe distributions 



measures the local slope of the xe distribution (NLD of an exponential distribution is constant). 
The NLD's extracted from the xe distributions of neutral pions associated with quarks and gluons 
as a function of xe is presented in Fig. 3. Because the xe distributions in the xe range above 0.8 
cannot be properly described by the fit (see Fig. 2), the NLD distributions are shown only below 
xe ~ 0.8. The NLD's for low xe regions (less than xe ~ 0.2) are much steeper than any other 
range, stay in lower values in the intermediate range (0.2 < xe < 0.7) and get larger for xe > 0.7 
for all quark and gluon trigger momenta. The NLD's from PYTHIA in the range of available 
PHENIX data are shown in the same Fig. 3 as open circles for each quark/gluon momentum. The 
mean xe extracted from the PHENIX data and the xe values used to extract "NLD" from PYTHIA 
by limiting the hadron acceptance of PHENIX are comparable to each other, as can be seen from 
Fig. 1 . In the latter comparisons of the PHENIX data to PYTHIA, the mean xe values labeled 
as "PHENIX acceptance" are used for higher momenta in PYTHIA where PHENIX data are not 
available. 

Finally, we have compared the exponential slopes of PHENIX data with PYTHIA where the 
used xe values are close to the mean values of PHENIX xe ranges in Fig. 4. They reveal the 
same increasing trend for the higher trigger pj t and the local slopes of quarks and gluons from 
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(2.2) 
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(a) quarks (b) gluons 

Figure 3: PYTHIA p + p at ^=200 GeV. Negative Logarithmic Derivative (NLD = ln(^)) of the 
xe distributions of 7r°'s resulting from quarks and gluons fragmentation in various quark/gluon momenta. 
As for the open circle labeled as "PHENIX xe range", the mean xe is calculated by limiting the associated 
charged hadron momenta from 1 to 4 GeV/c of PHENIX acceptance in the given trigger momentum. 
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Figure 4: Inverse xe slopes from PHENIX preliminary as a function of isolated y trigger momentum are 
compared with NLD's from xe distributions of K 0, s resulting from PYTHIA quark and gluon fragmentation 
with two set of xe ranges. The solid lines labeled as "PHENIX xe range", represent the local slopes in 
the expected xe value by considering PHENIX acceptance for charged hadron momentum in PYTHIA. The 
other shown as dotted line represents the NLD's for a fixed value of xe — 0.3 from PYTHIA. 
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PYTHIA are comparable to PHENIX data. The steepness of the exponential slopes in PHENIX 
data is driven by the experimental limit of the xe range, i.e away-side hadron can be measured 
only from 1 to 4 GeV/c in PHENIX, therefore xe value gets smaller as the trigger momenta gets 
higher. As for the intermediate xe (~ 0.3), the absolute slopes are slightly increasing as the trigger 
momentum gets larger and tend to converge into a constant above 10 GeV/c for quark and gluon 
fragmentation functions, shown as dotted lines. But the increment of the slopes for the gluon 
fragmentation functions (slope changes from 8 to 10) is rather larger about by factor of two than 
for the quark fragmentation functions (slope changes from 7 to 8) in the range 0.2 < xe < 0.7 (see 
Fig. 3) when the trigger momentum varies from 5 to 50 GeV/c. The fact that the absolute NLD 
values in all trigger bins agree with those from quark and gluon fragmenation function suggests 
that the measured xe distribution associated with direct photon is a representation of the quark or 
gluon fragmenation function. But a more detailed analysis with larger pj trigger bins and wider xe 
ranges is necessary to study trigger momentum and xe dependent effects which can influence the 
fragmentation function. 

3. Conclusions 

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured xe dis- 
tributions of charged hadrons associated with n° and isolated photon triggers in p + p collisions at 
■*/s = 200 GeV. It has been demonstrated in [6] that xe distributions measured in di-hadron correla- 
tions are insensitive to the fragmentation function. Direct photon associated spectra are expected to 
follow fragmentation function more closely, because in leading order two-to-two processes the trig- 
ger photon fixes exactly the momentum of outgoing parton by momentum conservation. However, 
several complications, like fragmentation photons and soft QCD radiation, are expected. Never- 
theless, the measured data presents a striking feature that slopes of the xe distributions, in both 
71° and photon trigger cases, continue to rise as a function of trigger momentum. We have stud- 
ied this rising trend in the case of photon trigger by fitting the same functional form as used in 
KKP parameterization [2] to the xe distributions of neutral pions resulting from fragmentation of 
u-quark or gluon in PYTHIA. As the fit function is not exponential, we use Negative Logarithmic 
Derivative (NLD) to measure local slopes. We found that the observed rising trend of slopes of xe 
distributions is driven by the experimental limits of the xe range. PHENIX can measure photons 
up to very large momenta while for hadrons the measurable range is more limited. Therefore the 
measured xe values for the higher trigger momenta are very small where the kj effect dominates 
and make the local slopes of the fragmentation function larger. Improvement of the statistical and 
systematic precision of the measurements with wider xe ranges and larger pj trigger should allow 
further tests of fragmentation function as well as understanding of kj effect quantitatively. 

References 

[1] T. Sjostrand, S. Mrenna, and P. Skands, JHEP0605 026 (2006). 

[2] B. A. Kniehl, G. Kramer, and B. Potter, Nucl. Phys. B597, 337 (2001), hep-ph/001 1 155. 

[3] A. M. Yu. Dokshitzer, V. Khoze and S. Troyan, Basics of Perturbative QCD (FrontiRres, 1991). 

[4] P. Abreu et al, Physics Letters B 398, 194 (1997). 



6 



Jet properties from direct y - hadron correlation in PHENIX at RHIC 



D.J Kim 



[5] R. D. Field and R. P. Feynman, Phys. Rev. D 15, 2590 (1977). 

[6] PHENIX, S. S. Adler et al, Phys. Rev. D D74, 072002 (2006), hep-ex/0605039. 

[7] F. Arleo, (2007), hep-ph/0701207. 

[8] I. Vitev and B.-W. Zhang, Phys. Lett. B669, 337 (2008), 0804.3805. 

[9] R. P. Feynman, R. D. Field, and G. C. Fox, Phys. Rev. D18, 3320 (1978). 

[10] J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 

[11] M. Nguyen, J. Phys. G35, 104121 (2008), 0805.1225. 



7 



